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Quantum anomalies, breakdown of classical symmetries by quantum effects, provide a sharp 
definition of symmetry protected topological phases. In particular, they can diagnose interaction 
effects on the non-interacting classification of fermionic symmetry protected topological phases. In 
this paper, we identify quantum anomalies in two kinds of (3-|-l)d fermionic symmetry protected 
topological phases: (i) topological insulators protected by CP (charge conjugation x reflection) 
and electromagnetic U(l) symmetries, and (ii) topological superconductors protected by reflection 
symmetry. For the first example, which is related to, by CPT-theorem, time-reversal symmetric 
topological insulators, we show that the CP-projected partition function of the surface theory is not 
invariant under large U(l) gauge transformations, but picks up an anomalous sign, signaling a Z 2 
topological classification. Similarly, for the second example, which is related to, by CPT-theorem, 
time-reversal symmetric topological superconductors, we discuss the invariance/non-invariance of 
the partition function of the surface theory, defined on the three-torus and its descendants gener¬ 
ated by the orientifold projection, under large diffeomorphisms (coordinate transformations). The 
connection to the collapse of the non-interacting classification by an integer (Z) to Zie, in the 
presence of interactions, is discussed. 
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I. INTRODUCTION 

Gapped states of quantum matter can be topologically 
distinguished by asking if given ground states can be adi- 
abatically connected to each other. Topologically equiv¬ 
alent phases can continuously be connected in the phase 
diagram (“theory space” in the language of the renor¬ 
malization group) with no gapless phase boundary nor 
quantum critical point separating them. On the other 
hand, topologically distinct phases are always separated 
by intervening gapless phases or quantum critical points. 

Quite often, it is meaningful to discuss theory space 
in the presence of a set of symmetries. Symmetries may 
prohibit the appearance of some phases distinct topo¬ 
logically from trivial phases, while they may also create 
a new topological distinction among quantum phases. 
Symmetry-protected topological (SPT) phases of mat¬ 
ter correspond to the latter case - they are topologically 
equivalent to trivial states of matter (such as an ionic 
insulator) in the absence of symmetries, while once a 
certain set of symmetries are imposed, they cannot be 
adiabatically connected to topologically trivial phases. 
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For partial references for recent works on symmetry pro¬ 
tected topological phases, see Refs. 1-11. 

As different SPT phases respect the same symmetry, 
they cannot be characterized by the Landau-Ginzburg- 
Wilson paradigm based on spontaneous symmetry break¬ 
ing. Hence one needs to look for an alternative organizing 
principle. Quantum anomalies, an intricate form of sym¬ 
metry breaking caused by quantum effects, have been 
proved to be useful in this context. Already in Laugh- 
lin’s gauge argument, a topological charge pumping pro¬ 
cess, i.e., the non-invariance of the system’s ground state 
under large U(l) gauge transformations, was used to es¬ 
tablish the stability of the quantum Hall states against 
interactions and disorder.For SPT phases, see recent 
works in Refs. 13-23. 

In this paper, basing on our previous works on (2-|-l)d 
SPT phases,the surface states of (3-|-l)d fermionic 
SPT phases are studied from the perspective of global 
quantum anomalies. We will discuss two examples: (i) 
the Dirac fermion surface state of (3-|-l)d bulk CP sym¬ 
metric topological insulators (TIs) and (ii) the Majorana 
fermion surface state of (3-|-l)d bulk reflection symmet¬ 
ric crystalline topological superconductors (TSCs). The 
bulk phase of the first example is a fermionic SPT phase 
protected by electromagnetic U(l) and CP [product of 
charge conjugation and mirror reflection (parity)] sym¬ 
metries. This example is CPT-conjugate to a (3-|-l)d 
time-reversal symmetric TI (class AH), and character¬ 
ized by a Z 2 topological number.^® The bulk phase of 
the second example is a fermionic SPT phase protected 
by fermion number parity and reflection (parity) sym¬ 
metry. It belongs to symmetry class D-I-R+ crystalline 
TSCs^® in Refs. 30-32. This example is CPT-conjugate 
to a (3-|-l)d time-reversal symmetric TSC (class DHI).^® 
At non-interacting level, class D-I-R+ crystalline TSCs 
in (3-|-l)d are characterized by an integer (Z) topologi¬ 
cal number (the Mirror Chern number), similar to their 
CPT partner, class Dili TSCs.^^ On the other hand, a 
number of recent works showed that the integral non¬ 
interacting classification of class Dili TSCs breaks down 
to Z 16 once interactions are included.Such col¬ 
lapses have been also reported in one and two spatial 
dimensions. 

Similar to our previous work,^® we enforce CP or re¬ 
flection symmetry on the surface theories by taking an 
orientifold projection.(See also recent discussion in 
Refs. 18 and 53.) In the first example, the resulting pro¬ 
jected theory is then shown to have global U(l) gauge 
anomaly. That is, the partition function of the pro¬ 
jected theory picks up a phase under large U(l) gauge 
transformations. This anomalous phase is shown to be 
a minus sign, and hence leads to the Z 2 classification. 
In the second example, by computing the global grav¬ 
itational anomaly®^’®® of the Majorana surface states 
of class D-I-R+ TSCs, we study the “collapse” of non¬ 
interacting classification. The resulting projected theo¬ 
ries are then shown to be anomalous under large diffeo- 
morphisms (coordinate transformations). 


In a similar vein, in Ref. 56, the (3-|-l)d Weyl fermion 
on the surface of the (4-|-l)d quantum Hall system is 
shown to fail to be modular invariant in the presence of 
a background U(l) gauge field. 

The rest of the paper is organized as follows. In the 
remaining part of this section, we introduce some no¬ 
tations that will be used in the main text. In Sec. H, 
we establish the gauge and diffeomorphism invariance of 
the (2-|-l)d Dirac fermion theory defined on a spacetime 
three-torus, following Refs. 57 and 58. In Sec. HI, we 
study (3-|-l)d TIs protected by electromagnetic U(l) and 
CP symmetries. The surface theory projected by sym¬ 
metries is shown to be anomalous, as its (projected) par¬ 
tition function is not invariant under large U(l) gauge 
transformations, but picks up a minus sign, character¬ 
izing the Z 2 classification of the bulk phase. (3-|-l)d 
TSCs protected by reflection symmetry are studied in 
Sec. IV, where we discuss the invariance/non-invariance 
of the surface partition function, defined on the three- 
torus and its descendants generated by the orientifold 
projection, under large diffeomorphisms. We then con¬ 
clude in Sec. V. 


Notations 


The partition functions of the (2-|-l)d surface theories 
discussed in the text can be represented in terms of par¬ 
tition functions of (l-l-l)d theories. Here, we summarize 
the properties of these (l-l-l)d partition functions. 

The partition function of a (l-l-l)d chiral fermion 
(Weyl fermion) on the two-torus with the modular 
parameter t = ti + iT2 G C, in the presence of spatial 
U(l) flux a and temporal U(l) flux b, is defined as^®’®® 






a-1/2 
6 - 1/2 


(0,t), 


( 1 ) 


where 77(r) is the Dedekind eta function and (u,t) 
is the theta function with characteristics. A^^j(t) has 
the following properties: 

Afa,b]i.T) = ^^+i.&](t) = 

<.](r + l) = e— 


The partition function of a (l-l-l)d massive Dirac 
fermion on with twisted boundary conditions (fluxes 
a and 6) is given by the “massive theta function” 

e[.,5](r;m):®o.®i 


0[a.h]('r;TO) 


_ g47rr2A(m;a) 


X n 


1 — e 


— 27rT2\/m^ + s^+27r2Tis+27rib 


(3) 
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where A(m;a) is the regularized zero-point energy: 


The Dirac fermion theory on three torus T® 


A(m;a) = x ^ 


sGZ+a 

OO pQO 


1 r 

2 I 


dk\/ m? + k^ 


-2^1: 


dte 


* cos(27rnaa;). (4) 




n—1 


The massive theta function 0[a,b](T; m) has the following 
properties: 


0[a,6]('r;TO) = 0[_a-&](r;m) = 0[a+r,b-hs]('r; m), 
0[a,6](T-f l;m) = 0[a.6+a](r;TO), 


0 [a.h] (-1 /t;w|t|) = 0 [b_a](T;m), 
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liin 0[a,&]('r;TO) = 

m—>0 




^fa.b]('r) 


r,sGZ, 

(5) 


II. LARGE U(l) GAUGE AND 
DIFFEOMORPHISM INVARIANCE OF (2+l)D 
FERMION THEORY 

In this section, we quantize the (2-|-l)d free Dirac 
fermion theory on a flat spacetime three-torus T^, in 
the presence of background U(l) gauge field and met¬ 
ric. The invariance of the partition function under large 
U(l) gauge transformations and 3d modular transforma¬ 
tions SL(3,Z), the mapping class group of T^, will be 
established. A discussion for the 2d modular invariance 
of the Dirac fermion theory on two torus as a warm 
up, is reviewed in Appendix A. 

We closely follow the analysis and notations in Ref. 
57. (See also Refs. 58 for related works.) In Ref. 57, the 
partition function of a chiral self-dual two-form gauge 
field on a 6 d spacetime torus T®, and its invariance un¬ 
der SL( 6 ,Z), the mapping class group of the six-torus, 
was studied. In Ref. 57, the theory is quantized (regular¬ 
ized) in a way manifestly symmetric under SL(5,Z). It 
was then shown that the partition function has an addi¬ 
tional SL(2, Z) invariance, and together with the SL(5, Z) 
invariance, the full SL( 6 ,Z) invariance was proven. By 
properly adopting this strategy, we show the SL(3, Z) in¬ 
variance and the large gauge invariance of the ( 2 -|-l)d 
Dirac fermion theory. 

While our focus in this section is on the complex or 
Dirac fermion, the case for real or Majorana fermions 
can be studied in a similar way. The modular properties 
studied here are expected to be straightforwardly gener¬ 
alized to higher dimensions, e.g., SL(n,Z) invariance for 
the partition function. 


1. Background metric 


A flat three-torus is parameterized by five “modular 
parameters”, Ri, 2 /Ro, a, j5, and 7 , where R^ are the 
radii for the /r-th directions, and a^j5,^ and related to 
the angles between directions 0 and 1 , 1 and 2 , and 0 and 
2, respectively. The dreibein is given by (^, A = 0,1, 2) 


/ Ro 0 OWl 0 0 \ 

0 i?i 0 -a 1 0 

\ 0 0 i?2 / V -7 -/3 1 / 

f Ro 0 0 \ 

I — (y.R\ R\ 0 1 , 

y -7R2 -PR2 R2 ) 


and its inverse is given by 


p* M _ 
— 


Ro 

0 

0 


Ro 

1 

Ri 

0 


0/3+7 

R2 


( 6 ) 


(7) 


such that and = 5^b- The Eu¬ 

clidean metric is 


A Be 
^ I'^AB 

/ Rl+ a^Rl + -aRj + p-fR^ -jRj \ 

= -aRl+pjRl Rl + P^Rl , 

y -7i?2 _p^2 ^2 j 

( 8 ) 


and the line element is given by ds^ = g^^dO^dO'', where 
0 < < 27r are angular variables. 

The group SL(3,Z) is generated by two modular 
transformations: 


Ui = 


/o 0 1 
10 0 
yo 1 0 


1 1 0 
1/2 = I 0 1 0 
0 0 1 


(9) 


The dreiben and metric are transformed as 
A (eL^)^^ = L/e^p, 

e*A {e*L-^)/ = 

ffppA {LgL^)^, = L/L,^gp,, (10) 

for any SL(3,Z) elements L = ■ ■ ■. In par¬ 

ticular. 


/ + 1)2^2 + (^ + ^3)2^2 1)3^2 + + 
iU 2 gU^)., = -(a - l)i?? + /3(7 + Rl + P^Rl -ml 

\ -{i + P)Rl -PRl Rl ; 


( 11 ) 
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which corresponds to the changes 

a —>-0 — 1 , 7—>-7 + ^ (while i?o, Ri, R 2 , and /3 are unchanged). 

The less trivial generator Ui can be further decomposed into two transformations as 


0-10 


1 0 0 


Ui = U[M, U[ = \ I 0 0 , M= 00-1 


( 12 ) 


(13) 


0 0 1 


0 1 0 


The transformation U[ acts on the metric as 


Rl + p^Rl aRl - p-fRl 
-iRl 


I . 11-1 -TM -‘'-2 PR 2 

{U[gU[^)^^= aRl-MRl Rl + a^Rl + i^Rl -jRl 


PRl 


Rl 


(14) 


which corresponds to the changes 
i?0 —>■ Ro/\'r2d\, Rl Rl\T2d\, O' 

with 

while M acts on the metric as 

M 


-a/\T 2 df (or T 2 d-t-l/r 2 d), 7 -/3, -t 7 (while i ?2 is unchanged), 

(15) 


T2d = o + ir 01, = R^/Rv, 


(M 5 M% = 


i?Q + a^R\ + 7^i?2 iRl —oR\ + (i^R 

iRl Rl PRl 

-aRl + P"fRl ml Rl + P'^Rl 


(16) 


(17) 


The Euclidean action for the Dirac fermion on is then given by 
1 


Sr — 


(2^)2 

1 




^27ri?o p2'KR\ 

dr / 


/•27ri?2 


dx 


dyip 


+ (a/3 + i)^T°dy + T^S, + + T^dy 

Ko ^ 


Rn 


i?i 


pJ, 

(18) 


where ip is the two-component Dirac field, ip = ip'^T^, where {ar,ax,ay) = a represents the background U(l) 

T = i?o 0 °, X = Rm, y = R 20 ^, and the gamma matrices gauge field twisting the boundary conditions. 

T"^ satisfy {r^,r®} = 25^^. 


2 . Background flux 

In addition to the background metric, we also intro¬ 
duce the background U(l) gauge field (flux) on to 
twist the boundary conditions of the Dirac fermion the¬ 
ory. More specifically, in the path integral language, we 
consider the boundary conditions 

ipir, X + 27ri?i, y) = x, y), 

ipir, x,y + 2nR2) = x, y), 

'ip{T + 27ri?0i X — 27rai?i, y — 27r(a/3 -I- 7)i?2) 

= e2-“^^(T,a;,y), (19) 


3 . Partition function 


We now quantize the (2-|-l)d theory and compute 
the properly regularized partition function, denoted by 
Z[a] ( 5 ) j which depends on the background flux a and 
metric g. The partition function can be evaluated by 
the path integral on T^, .^[a](y) = /,'(/)] exp(—S'£;), 
with ip satisfying the twisted boundary conditions (19), 
or alternatively, in the operator language, by the trace 


Zl3]ig) = Tra^a, 


^2Tri(a^-l/2)F 


^-2-kRoH' 


( 20 ) 


where H' is the ’’boosted” Hamiltonian (in the presence 
of non-vanishing angles a, /3, and 7 ) obtained from Se 
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and given by 

H'= H-ia^P,-i{aP + -i)^Py, (21) 

iXO riQ 

with 

H = J dxdy^ (r^d, + Id^r^dy + V', 

Pi = J dxdyip^{-iditp), i = x,y, (22) 

being the Hamiltonian and momenta. Tia^^ay means the 
trace is taken over the Fock space of the fermion theory 
for the spatial boundary conditions specified by Qx and 
tty. The twisted boundary condition in the r-direction 
is implemented by an operator insertion exp[27ri(ar — 
1/2)^], where F is the fermion number operator. 

The fermion held operator satishes the canonical anti¬ 
commutation relation 

{V'a(r)>^^(''')} = (27r)^da/3 

mi ,m 2 GZ 

X 5{x — x' + 2'KmiRi)5{y — y' + 27rm2i?2)) (23) 


where r = {x, y) and a, /3 are spinor indices. The trace 
can be evaluated explicitly by the Fourier mode expan¬ 
sion of the fermion held operator. With the twisted 
boundary conditions, the fermion held operator is ex¬ 
panded as 


V'(r) 


1 

\/ i?ii ?2 


E E 

Sy^Wl-\-CLy 


(24) 


where s = (s^, Sy) and 

{^a(s), V'^(s')} = darks'- (25) 


Correspondingly, the Hamiltonian can be expanded as 


i?= E E ^^(s)'H(s)v;(s), 


H(s) = r° 



Rl i ?2 i ?2 


(26) 


The single-particle Hamiltonian 'H(s) can be diagonal¬ 
ized with eigenvectors u±{s) and eigenvalues ±e(s): 


The Hamiltonian can be diagonalized by the eigen basis 
x(s) := [x-i-(s),x_(s)]^, which are related to the original 
fermion operators ^(s) as 


^l(s) 


Wl+(s) Ui_(s) 

X-h(s) 

>2(s) . 


_U2+{s) U2_(s)_ 

.X-(s) . 

X-k(s) 


ul+{s) u5+(s) 

^l(s) 



_m):_(s) u*_{s) _ 

.^2(s) . 


(29) 


The Hamiltonian in the eigen basis is given by 
H = ^e{s) x+(s)X-h(s)- xUs)x_(s) 

S 

= E^(^) [xUs)X-h(s)+X-( a;)xUs)J -E^(s) 

S 5 

=: H : +Egs (30) 


where : • • • : is the normal ordering with respect to the 
Fock vacuum and Eqs = ^he ground-state 

energy. 

The ground-state energy needs to be properly regular¬ 
ized. As shown in Appendix B, we have 


A, 


'GS[; 


= "E 


92S^Sj 


= ^Vdet(52z,) E 


cos ( 27 rio‘nij 


n^ 0 GZ 2 




(31) 


where Sx^y are separated into their integral and fractional 
parts as 


Si = ni + ai, rii G Z, i = x,y. (32) 


Similarly, the ground-state momentum and the fermion 
number can be regularized as 

(-P^GS = E ^ E + E 

Si Si>0 Si <0 

= C(-l,a,)-C(-l,l-aO = 0, 

Fgs=Ei-Ei = 0 ’ (33) 

5 S 

where C('Si 2 ;) = + Hurwitz zeta func¬ 

tion dehned by analytic continuation from the region 
Re(s) > 1 and we have used C{—l,x) = ^ — ^[x—^) . 

With the regularization, the partition function (with 
boundary conditions twisted by Ux and Uy), given by the 
trace in (20), is evaluated as 


'H{s)u±{s) = ±e{s)u±{s), 

£(S) = 


Z[a]ig) = e 


_ „ —27r/JoEGS 


n n 

Sy^'Z~\~CLy 


1 — e 


— 27rRQe{s)-\-27riaSj.-\-27ri{a(3-\-'y)sy-\-27riaT 


(34) 


where 


in 

92 


911 

921 


9l2 

922 



1 



Rl 


/3 

+ 


Rl 


which can also be expressed as the infinite product of 
massive theta functions defined in (3): 

•^[a](5)= ^[a^+Psy,ax+^Sy]{.'^2d'-!'ri2Sy) ■ (35) 

Sy 


(28) 
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We now demonstrate the invariance of the partition 
function under large U(l) gauge transformations and 
modular transformations. 

Jf. Large U(l) gauge invariance of the partition function 

We first check the invariance of the partition func¬ 
tion under large U(l) gauge transformations 
0‘x,y,T + 1- The invariance under ax,T —t tix,T -I- 1 is ob¬ 
vious from Eq. (35), using the properties of the massive 
theta function listed in (5). To check the invariance of 
the partition function under ay ^ ay + 1, we note that 
this amounts to a simple shift Sy —> Sy -I- 1 in the infinite 
product in Eq. (35). To sum up, we conclude the large 
U(l) gauge invariance of the partition function. 

5. Modular invariance of the partition function 

By using the expressions (34) and (35) of the parti¬ 
tion function, we can show that (g) has the following 
property: 

Z[LB\{LgL^) = ^[a](5)> 
or Z[a](LgL'^) = Z[i-ia](g), (36) 

where L S SL(3, Z). This means the Dirac fermion, when 
coupled to both background U(l) gauge field and met¬ 
ric, is anomaly-free under any large diffeomorphisms (to¬ 
gether with the induced gauge transformations) on T^. 
The claim (36) can be shown by checking how Z[a](g) 
transforms under Ui = U[M and U 2 , defined in Eqs. (9) 
and (13). Here we leave the detail of the derivation to 
Appendix C. 

We now show that the partition function, once pro¬ 
jected by the fermion number parity only [i.e., in the 
absence of U(l) gauge field], is modular invariant. We 
consider the sum over all periodic/antiperiodic boundary 
conditions, which are specified by boundary conditions 
a = ( 0 , 0 , 0 ), ( 1 , 0 , 0 ), ..., ( 1 , 1 , 1 ), corresponding to the 
23 = 8 spin structures when defining fermions (spinors) 
on . The resulting total partition function is given by 

^‘°‘(5)= E (37) 

Ox,y,x— 0,1/2 

where Cg are weights (“discrete torsion”) assigned to dif¬ 
ferent sectors with partition functions twisted by a. Erom 
Eq. (36), we see that by choosing Cg = 1 for all a, the total 
partition function is modular invariant: 

Z^°*{LgL^) = Z*°*{g), L G SL(3,Z). (38) 

III. SURFACE THEORY OF (3+l)D CP 

SYMMETRIC TOPOLOGICAL INSULATORS 

Based on the result from the previous section, now we 
can compute quantum anomalies of an anomalous surface 


theory and interpret them as a signal of the existence of 
the nontrivial bulk SPT phases. In this section, we iden¬ 
tify a global U(l) gauge anomaly of the surface theory 
of (3-|-l)d CP (charge conjugation x reflection) symmet¬ 
ric TIs, which are related to, by CPT-theorem, (3-|-l)d 
time-reversal symmetric TIs. 

A. Surface theory 

Let us consider the following surface Dirac Hamilto¬ 
nian of (3-|-l)d TIs: 

^^(2^ y c^a:(i2/V'^(r)(-zCT292, - *cri5y)'0(r)> (39) 

where cri. 2,3 are the Pauli matrices, the spatial coor¬ 
dinate r = {x,y) G [0,27ri?i) x [0,27rii2) parameter¬ 
izes the 2 d surface, and tf{r) = [ijji{r),ip 2 {t)]'^ and 
' 0 ^('') = ['*/’!('')I' 02 ('')] are two-component fermion annihi¬ 
lation and creation operators, respectively. The Hamilto¬ 
nian is invariant under the following time-reversal sym¬ 
metry: 

= *CT 2 V'(r), = (-1)^, (40) 

where F is the fermion number operator. This time- 
reversal symmetry prohibits the mass term since 

Alternatively, in the following, we will take the Dirac 
Hamiltonian (39) as a surface theory of bulk CP sym¬ 
metric TIs. By CPT-theorem, the classification of CP 
symmetric TIs are expected to be the same as the clas¬ 
sification of time-reversal symmetric TIs. That is, CP 
symmetric insulators in (3-|-l)d are also classified by a 
Z 2 topological number.^® Within the surface theory, the 
action of CP symmetry is given by 

= a3[tfj\x,2TTR2 - y)f, 
{x,y){‘l^^)~'^ = i{j{x,27rR2 - 
= 1. (41) 

[This is the only CP symmetry of the Dirac kinetic term 
'H(kx,ky) = kx(T 2 + kyai since {—kx,ky)a^^ = 

—^.{kxyky).] Fermion bilinears in the surface the¬ 
ory are transformed as M'if —>■ = 

—'01In particular, the mass is odd under CP; 
the surface theory, at least at quadratic level, cannot be 
gapped without breaking symmetries, U(l) x CP. On 
the other hand, a CP preserving mass exists if we dou¬ 
ble this theory (or more generally if the number of the 
Dirac fermions is even), and the corresponding surface 
theory can be gapped. Since massive fermions can al¬ 
ways be regularized to construct a well-defined quantum 
theory, an even number of the surface fermions (39) is 
always anomaly-free (while preserving the symmetries). 
However, an odd number of the surface fermions may suf¬ 
fer from anomalies. In the following, we will identify a 
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quantum anomaly of the surface theory (39) under large 
U(l) gauge transformation when CP symmetry is strictly 
enforced. 


B. Projected partition function by CP symmetry 


We now consider CP projection of the surface theory 
(39) and ask if the projected theory is still invariant un¬ 
der large gauge transformations (as we have seen in the 
theory without symmetry projection). This leads to for¬ 
mulating the fermion theory on unorientable spacetime 
manifolds such as 5”^ x K, where K is the Klein bottle. 
As our main focus here is on large gauge transformations 
but not on modular transformations, the parameters a, 
/3, and 7 are set to zero in the following discussion. Also, 
the twisted boundary conditions by U(l) are consistent 
with CP symmetry only when ar^x = 0,1/2, while ay is 
not constrained by CP symmetry. We will thus study the 
large gauge transformation —>■ Oy -I- 1. 

The partition function of the surface theory projected 
by CP is given by the trace 


Tr, 


Ora; 


1(1 _|_ <^_^)g27rhaT-l/2)Fg-27rfloff 


(42) 


Upon projection by CP, we will focus on the CP sym¬ 
metric boundary conditions Ot-.x = 0,1/2 and ay G [0,1). 
The first term in the projected trace, which is invari¬ 
ant under Oy —>■ Oy -I- 1, is already discussed in Sec. II 4. 
Our focus below will be the second term in the projected 
trace, which we call the CP twisted partition function: 

ZgP = Tia^^ay ^. (43) 

Because of CP symmetry, from the eigenvectors u± at 
s, we can construct eigenvectors at s = {—Sx,Sy): 

■H(s)(T3?4(s) = =F£(s)cr3ui(s)- (44) 


The CP operator acts on the Fourier components of the 
original fermion operators as 




0-3 



(45) 


On the other hand, the CP action on the eigen basis 
X± , x± is deduced as 

_ (u+(s)|cr3A:|w+(s)) (M+(s)|cr 3 A:|M-(s)) 1 r t, .tT 

(u_(s)icr3A:iM+(s)) (M_(s)icr3A:iu_(s)) J i ''-I ’ 

(46) 

where K is the complex conjugation operator, and 
(M±(s)|cr3Ar|M±(s)) = u^{s) ■ tT3?I5_(s), etc. Since u±{s) 
and (T3'it5-(s) are both eigenvectors of Hamiltonian 77(5) 
but with different energies, their overlap should be zero. 
Therefore, 

{^^)x+is){^^)~^ = {u+{s)\a3K\u-{s))xiis), 
{'^3^)X-{s){^ = (M_(s)|cr3A:|u+(s))xl(s). (47) 


The transformation law of x± under CP depends on the 
choice of eigen functions u±. A choice for the eigenvectors 
is 


i.as)) ^ /j 


±AV|A| 

1 


A = 


(48) 


For this choice of eigenfunctions, 

(u+(s)|cr3A:|M_(s)) = (■u_(s)|ct 3 A:|m+(s)) = -1. (49) 

As we can choose the phase of the eigenvectors freely, 


1 




is also an eigenfunction. In this gauge, 

(w+(s)|cr3A:|M_(s)) = (u_(s)|cr3A:|M+(s)) = 1. 


(50) 


(51) 


In either choice, the result can be summarized as 

Crf ^)x+{s){'rf = V+X-i^): 

y yx-^y = 11-X+y, (52) 


where r]± is s-independent. The product rj := rj+ri- = 1 
is gauge invariant. 

The CP twisted partition function can then be 
computed explicitly as 


•CP 

[a 


— 2-kRqEgs P, , 

[aa:,ay\ 


X n n 


1 


g-4ir/Joe(s) 


(53) 


where the prefactor Pia^^ay] is the CP eigenvalue of the 
ground state (the Fock vacuum). Note that the partition 
function does not depend on Ot, which is projected out 
by CP.®^ In the following, we denote = Z^^ and 
consider the two cases Ux = 0,1/2 separately. 

a. Periodic boundary condition in the x direction, 
Ox = 0; In this case, we can factor the twisted parti¬ 
tion function into the product of 2d massless (s^ = 0) 
and massive {sx ^ 0) modes 


^CP 


'[ay:=0,ay] 

27ri?o^. £:(Sa:;=0,ST, 

X e ^ 


(54) 


JJ" ^ _ g-4ir/Joe(sj,=0,Sa) 

Sy ^'Z-\-CLy 


X n 

sxez+ 


4TrRo J2s„ =(=) 




Sy 


which can be expressed (as the sum J2sy ^(^) i® regular¬ 
ized) in terms of (l+l)d partition functions defined in 
Sec. I: 


yCP _ p _ 

'^[aj,=0,a„] “ 7“a:=0,a»]® 


— 7rz(ay —1/2) 


X n 0K.o](2fro2;r2is^), (55) 

Sx 




















where = R^/R^. When imposing CP symmetry, it is 

reasonable to assume the CP-eigenvalue does not change 
under ay ^ ay + i.e., P[a,,=o.a„+i] = P[a^=o,ay]- Un¬ 

der this assumption, we have 

7CP —_yCP 

^[a^=0,ay + l] ^[a^=0,ay]^ 

where we note both and 0[a„,o] invariant under 

fla; —t Oa: + 1- The anomalous minus sign under the large 
gauge transformation, which comes from the 2d massless 
modes {s^ = 0) but not the massive modes (s^ ^ 0), 
signals a Z2 topological classification: the CP projected 
theory can only be realized as the surface theory of a 
(3-|-l)d bulk CP symmetric TI, which is CPT-conjugate 
to a (3-|-l)d time-reversal symmetric TI.^® 

b. Antiperiodic boundary condition in the x direction, 
ax = 1/2.' In this case, the twisted partition function is 
given by 

^[ar=l/2,a„] = U[a,, = l/ 2 .a„] H ©K .0] (2ir02 ; ?'21 Sa;). 

s,ez+-i/2 

(57) 


Observe that there are no 2d massless modes arising in 
the expression of ^[^^^1/2 ay] (while the product of all 
massive modes changes from Y{sy&+ to ns„GZ+-i/2)- 
This partition function is anomaly-free under the large 
gauge transformation, i.e.. 


7CP _ yCP 

'^K = l/ 2 .a„-|-l] — ^[a,,=l/ 2 ,a„]' 


(58) 


For arbitrary number N of the Dirac fermion flavors, 
the result is summarized as: 

. (59) 

The surface theory, as projected (or twisted) by CP, is 
anomaly-free if and only if = 0 mod 2. This charac¬ 
terizes the Z2 classification of the bulk SPT phase. 


IV. SURFACE THEORY OF (3-|-l)D 
REFLECTION SYMMETRIC CRYSTALLINE 
TOPOLOGICAL SUPERCONDUCTORS 


In this section, we identify a global gravitational 
anomaly of the surface theory of (3-|-l)d reflection sym¬ 
metric crystalline TSCs, which are related to, by CPT- 
theorem, (3-|-l)d time-reversal symmetric TSCs. While 
the Z2-type (gauge) anomaly in the surface of CP TIs 
agrees with the non-interacting classification of the bulk 
phase, the (gravitational) anomaly in the surface of re¬ 
flection symmetric TSCs, as we will discuss later, sees 
only the reduction of non-interacting classification, and 
hence can detect the effect of interactions (in the case 
that the bulk gap is not destroyed by the interactions). 


A. Surface theory 

At the quadratic level, time-reversal symmetric super¬ 
conductors in symmetry class Dill are classified by an 
integer topological invariant, the 3d winding number vA 
The topological invariant counts the number of gapless 
surface Majorana cones. For example, the B-phase of 
®IIe is a TSC (superfluid) with v = 1, and hosts, when 
terminated by a surface, a surface Majorana cone, which 
can be modeled, at low energies, by the Hamiltonian 

H = J d‘^rX^{-ia3dx-i(7idy)X, (60) 

where cri_2,3 are the Pauli matrices, the spatial coordi¬ 
nate r = {x,y) G [0, 27ri?i) x [0, 27ri?2) parameterizes the 
2d surface, and A(r) is a two-component real fermionic 
field satisfying A^(r) = A(r). The surface Hamiltonian is 
invariant under time-reversal tA' defined by = 

icr2A(r), where is equal to the fermion number par¬ 
ity = (— 1 )'^ = j2^ / '^^rA^tT2A. For TSCs with 
V = Nf, the surface modes can be modeled by iV/ copies 
of the above surface Hamiltonian. 

While, at the quadratic level, , one can verify that, for 
an arbitrary integer n = Nj, surface Majorana cones 
are stable against perturbations the surface Majorana 
cones may be destabilized once interactions are included. 
A number of arguments, such as “vortex condensation 
approach”, “symmetry-preserving surface topological or¬ 
der” , “cobordism approach” and so show that 

the surface Majorana cones are unstable against interac¬ 
tions when 1 ^ = 0 mod 16, reducing the non-interacting 
integer classification to Zie. 

Here, instead of time-reversal symmetry, we consider 
its CPT-conjugate, reflection or parity symmetry, which 
acts on the Majorana field as 

= a^\{x, 2 'KR 2 — y), = (61) 

Upon demanding the invariance under parity (61), the 
Majorana Hamiltonian (60) describes the surface of sym¬ 
metry class D -I- R+ crystalline TSCs, which are, at the 
quadratic level, classified by the integral mirror Chern 
number.Based on CPT-theorem, we expect, upon 
the inclusion of interactions, the integer classification col¬ 
lapses down to Zig. 

To see the stability of the gapless Majorana mode at 
the quadratic level, note that the mass A^(T2A is odd 
under parity (61) and prohibited. It is also interesting 
to note that while the uniform mass is not allowed, one 
could consider f m(r)A^cr2A with m{x, 2 'KR 2 — y) = 
—m{x,y). This perturbation gaps out the most part of 
the surface, but not completely. At the fixed points of P 
symmetry, y = 0 and y = 'XR 2 , it leaves gapless modes 
localized at the domain walls. Note that this is similar to 
the chiral mode localized at a mass domain wall on the 
surface of time-reversal symmetric TIs. The difference, 
however, is that in the present case, the mass domain 
wall, as a whole, preserves the reflection symmetry, while 
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the domain wall on the surface of TIs breaks time-reversal 
symmetry, except at the domain wall. The gapless mode 
at the domain wall consists of the Nf copies of Majorana 
fermions propagating in either +x or —x directions, de¬ 
pending on the overall sign of the mass domain wall for 
each flavor. For even Nf, we can always choose (tech¬ 
nically) a set of mass parameters such that the gapless 
modes at the domain wall are made nonchiral [e.g., by 
choosing different signs of the masses for different flavors 
of Majorana fermions). In this case, reflection symmetry 
acts on the (l-l-l)d gapless domain-wall fermions as an 
unitary on-site Z 2 symmetry. Using the result in Ref. 24, 
it can be shown that such gapless domain-wall states can 
be gapped without breaking the symmetry if the number 
of the non-chiral states is 0 mod 8. This means, when 
Nf = 0 mod 16, we can gapped out the surface of the 
crystalline TSCs while preserving the reflection symme¬ 
try at the same time. This gives the Zig classification, as 
expected to the same as the case of class Dill TCSs, of 
the class D -I- R-l- crystalline TSCs, upon the inclusion 
of interactions. A similar argument for the Zg classifica¬ 
tion of interacting crystalline TIs protected by reflection 
symmetry can be found in Ref. 64. 


B. Projected partition function by 
refiection/parity symmetry 


We now study the presence/absence of (global) gravi¬ 
tational anomalies of the surface theory (60), which sig¬ 
nals the existence of the nontrivial bulk SPT phases. For 
convenience, we double the degrees of freedom and con¬ 
sider Dirac instead of Majorana fermion fields. (This is 
purely a matter of convenience. The analysis below can 
be repeated without referring to the Dirac fermion, and 
can be done solely in terms of the Majorana fermion.) 
The number of Dirac fermion flavors will be denoted by 
N, which corresponds to Nf = 2N in term of the original 
Majorana fermions. 

Our starting point is the partition function = 

with H' given by Eq. (21). 
Let us now include the effects of parity symmetry by 
including twisted boundary conditions by parity. First, 
note that the modular parameters /? and 7 are odd under 
parity. Hence they will be set to zero henceforth, fi = 
7 = 0, to consider the parity twisted partition function. 
While SL(3, Z) acts on the metric g = g{Ri, a, /3,7), there 
is an SL(2,Z) subgroup generated by U[ and U 2 , acting 
on the “reduced” set of the modular parameters, gp = 
gp{Ri,a) = g{Ri,a,/3 = 7 = 0).®^ With the reduced 
set of modular parameters by parity, the total partition 
function, which is generated by projection by parity and 
the fermion number parity, is given by 


E ^G[ZG{9p)f, 

GeSG3 


ZG{gp) = Trg^ Gj, 




(62) 


(0 T (it) 




FIG. 1. The three-torus and its (unorientable) descen¬ 
dants generated by the orientifold projection. While the y- 
boundary condition is twisted by Gy = 1 or 1#/, the r- and x- 
boundary conditions are twisted by {Gt,Gx) = , 1#^““'), 

(_^^2a, ,^2a,,), ^ ^ ^ ), aS shown 

in figures respectively. (Un)twisted boundary condi¬ 

tions are represented by arrows with the same color. 


where SG = ^^^f} is the symmetry group of 

the surface fermion theory, and ec are weights assigned 
to different sectors with partition functions twisted by 
G = (Gt, Gx, Gy), where for each direction, the boundary 
condition is twisted by G^ = \,'^f,SZ, SZtSf. 

Not all sectors of the total partition function are 
mixed by SL(2,Z). We can then divide different sectors 
into groups, and study the action of SL(2,Z) on each 
group separately. In the following, we will focus on the 
sectors generated by twisting y-boundary condition by 
t,tSf, and by twisting r- and x- boundary conditions by 
\,tSf,SZ,tSf^. For a given j/-boundary condition, there 
are 4^ = 16 sectors in total, and the corresponding par¬ 
tition function is 

4d(5P)= E (63) 

GgSG2 

where G = (Gx,Gx), and Qy = 0,1/2 represents the 
y-boundary condition. We will consider the cases of 
Oy = 0 and Gy = 1/2 separately, as they are not 
mixed by SL(2,Z). The remaining sectors can be gener¬ 
ated by twisting j/-boundary condition by 3^ and ^t^f. 
Twisting by these group elements gives rise to what 
can be interpreted as ’’open sectors” (partition func¬ 
tions on orbifolds) as noted by Horava.^® In this pa¬ 
per, however, we will focus on the 32 “closed” sec¬ 
tors generated by twisting with Gy = 1,1^/. The re¬ 
sulting closed orientable/unorientable three-manifolds, 
where the (twisted) partition functions are evaluated, are 
shown in Fig. 1. 

In the following, we present the analysis of the twisted 
partition functions for the case of Oy = 0. The detail 
of the calculations is left to Appendix D. The analysis 
for the case of Oy = 1/2 is similar and in fact simpler. 
In short, for Oy = 1/2, the total partition function (63) 
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can be made modular invariant for any number of Dirac 
fermion flavors, N. See Appendix F. 

On the other hand, the total partition function for 
tty = 0 can or cannot be made modular invariant, de¬ 


pending on N. For Uy = 0, there are 16 sectors in 
total, generated by twisting by and in the t 

or/and x directions. We divide these 16 sectors into four 
sets (t-tv) by {Gr,G,) = “=), 

(^ 2 a^, ^ ^ ^ respectively. (There 

are four sectors in each set.) The symmetry-twisted par¬ 
tition functions ^ ^ 2 ay ^ for each set are then given 

by (see Appendix D) 

X[a^,a^]i9p) = (64) 

X[l^,a^]i9p) = ^^,,a,]{'^2d)-4^^_a^_i](T-2d)0[^^_2a,]> 

X[Z,a^]{9p) = ^K,a,]('^ 2 <i)A^^_i_^^j(T 2 d) 0 [ 2 a^_„^], 
X[a^,a^]i9p) = ^^,,,a,]('^ 2 <i)A^^_i_^^_lj(T 2 d) 0 [ 2 a^_„^_„^], 

where ar,x = 0 , 1/2 and we have introduced the functions 
(^ 2 d;ri 2 ) by 

0L „ , 


0 / 


0 


tit 


0 f: 


n [ 0 K,a,] (T- 2 d;ri 2 Sy)] , 

,a-r] (2'?'2(i5 ^ 12 'Sy); 

,ax] ('^ 2 (i/ 2 , 2 ri 2 Sy), 
,ax]('^2d/2 + l/2;2ri2Sy). 


n Q 

s„GZ+ 

n 0 

s„ez+ 

n Q 

SyGZ+ 


(65) 


When evaluating the partition sum (63), constant pref¬ 
actors may show up, but are not displayed in the expres¬ 
sions (64). These prefactors correspond to parity eigen¬ 
values of the ground states in different sectors [which 
might depend on the modular parameters and fluxes but 
are assumed to be SL(2,Z) invariant], and can be ab¬ 
sorbed to the (redefined) weights eG,a„ in Eq. (63). 

We now ask, for a specific choice of N, by sum¬ 
ming these partition functions with some set of weights, 
if we can construct a modular invariant. The trans¬ 
formation properties of the twisted partition functions 
X[a™a ] under SL(2,Z) (generated by U[ and U 2 ) can 
be deduced from the properties of and 0 shown 

in Sec. I; see Appendix E. It can be shown that if 
and only if IV = 4n (n = 1,2,3,...), i.e., Nj = 8n, 
a modular invariant can be constructed. In addition, 
while SL(2,Z) invariance can be achieved for N = An, 
there is a distinction between n = 2k — 1 and n = 2k 
{k = 1,2,3,...), i.e., N = 8k - A {Nf = 16k - 8 ) and 
N = 8k {Nf = 16k). To be explicit, the twisted par¬ 
tition functions in set (i) are closed under SL(2,Z) and 
a modular invariant can be constructed for any N. For 
the twisted partition functions in set (ii — iv), we con¬ 
sider a weighted sum where 


Xt = Xi = xfi,o]> Xi = Xp,i], X4 = Xjl,i]- 

When N = 8k — A, the SL(2,Z) invariance is achieved 
when 




1 5 *=2 5 ’ ^4 ’ ^1 ’ ^2 5^3 ’ ^^4 ? 5 *=2 ’ ^3 ’ ^4 


) ( 66 ) 


— (Olj U2 j 03j 0,3, Oi, 03, 02, 03, —Oi, —03, —03, —02), 


where ai=i_ 2,3 are arbitrary phases (signs). Thus, when 
N = 8k —A, the trivial choice, ef = l for all {A, i), is not 
allowed. When N = 8k, on the other hand, the SL(2,Z) 
invariance is achieved when 


(e 


a it it ii Hi Hi Hi Hi iv iv iv Hv\ 
1 5 ^2 5 €3 5 £4 ? 5^2 5 £3 , €4 , , 62 , €3 ,64 j 

= (oi, 02,03, 03, Oi, 03,02,03, Oi, 03,03, 02). 


(67) 


Dimensional reduction to the edge theory of the (2+l)d 
fermionic SPT phase with Z 2 x Z 2 symmetry. The 
SL(2,Z) invariance for N = An {Nf = 8n) may be 
understood by taking the limit i ?2 —t 0 (ri 2 —?> 00 ). 
In this limit, all massive theta functions become 1 and 
the total partition function constructed here reduces to 
the form of the (l-l-l)d partition function projected by 
Z 2 X Z 2 symmetry (fermion number parity conservation 
for each chirality), which is the edge theory of the ( 2 -|-l)d 
SPT phase with spin parity conservation.^^ In the latter 
case, the SL(2,Z) invariance of the Nf = 8n symmetry- 
projected partition function indicates that 8n helical Ma- 
jorana edge modes [in ( 1 + 1 ) dimensions] can be gapped 
without breaking the Z 2 x Z 2 symmetry. 


V. DISCUSSION 

We have studied global anomalies on surface theories 
of (3+l)d topological insulators and superconductors. 
For CP symmetric TIs, which are related to, by CPT- 
theorem, time-reversal symmetric TIs, there is a global 
U(l) gauge anomaly if the number of the surface Dirac 
fermion is odd, characterizing the Z 2 classification of the 
bulk phase. For reflection symmetric TSCs, which are 
related to, by CPT-theorem, class Dill TSCs, a global 
gravitational anomaly is present in the surface theory 
when Nf 7 ^ 0 mod 8 . The corresponding bulk state is 
topologically distinct from trivial states of matter even in 
the presence of interactions, as far as the bulk gap is not 
destroyed by the interactions. On the other hand, the 
weights cq, determining the relative weights among par¬ 
tition functions in different sectors, have 16-periodicity 
as a function oi Nf. Our analysis thus presents an al¬ 
ternative approach to the collapse of the non-interacting 
classification. 

For the cases where we do not find any inconsistency 
(anomaly), i.e., the case of TSCs with Nf = 8, the sit¬ 
uation may be more subtle. First of all, the theory may 
suffer from other forms of inconsistency, which have not 
been studied here, and hence particular calculations pre¬ 
sented in this work does not immediately conclude that 
the corresponding (3+l)d bulk theories are topologically 
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trivial. Recall that we have not included the partition 
functions twisted in the y direction by ^ and [see 

comments below Eq. (63)]. 

Moreover, we studied the problem of global anoma¬ 
lies by considering surface theories on (and its de¬ 
scendants generated by the orientifold projection). Even 
when the theory is shown to be consistent on it may 
be anomalous once put on a different three-manifold. The 
situation is better understood for 2 d conformal field the¬ 
ories (CFTs), where once the consistency of the theories 
at genus one (torus) is established, they can be consis¬ 
tently defined on any (oriented) Riemann surfaces. Eor 
3d CETs, there is no such known fact. Eor this reason, 
our quest for anomalies in the surface theories may not 
be complete. Nevertheless, our study on anomalies of 3d 
massless fermions has shown some interesting and non¬ 
trivial results.®® 

Finally, it is interesting whether our approach can 
be related to the gapped surface states of (3-|-l)d SPT 
phases that develop symmetry-respecting topological or¬ 
ders. Such connection is recently investigated in Ref. 23 
in the case of the SU(2) global anomaly.®^ Extending such 
connection to a generic set of interacting SPT phases is 
left for future studies. 
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Appendix A: The Dirac fermion theory on two-torus 


In this appendix, we review the modular invariance, 
the SL(2,Z) invariance, of the Dirac fermion theory on 
two-torus T^. 

For a flat T^, the zweibein can be factorized as 


e 


A 




i ?0 0 

0 i?i 


1 0 
—a 1 


Rq 0 
— (y.R\ R\ 


, (Al) 


and its inverse is given by 


(A2) 

such that = 6^'' and = 5^b- Here i?o 

and Ri are the radii for the directions 0 and 1 , and a is 
related to the angle between the directions 0 and 1. The 
Euclidean metric is then given by 


4^ = 


1 a 

Ro Ro 

0 ir 


= 


^ i 


Sab 


Rl + a'^Rl -aRl \ 
-aRl Rl ) ’ 


(A3) 


and the corresponding line element is 

ds^ = gf^^dO^dO^ = Rl{de°y + R^idO^ - ad6°f, (A4) 

where 0 < < 27r are angular variables. 

The group SL(2,Z) is generated by two transforma¬ 
tions: 


t/i = 


0 -1 
1 0 


[/,= 


1 1 
1 0 


(A5) 


SL(2,Z) transformations on the zweibein and metric are 
induced by 


e^„, A (eL^r ^ = L 


r j ^ — L/e^p, 




* M 




iLgL^)^^ = LpPL,^gp,, (A 6 ) 

for any SL(2,Z) element L = • • •. In partic¬ 

ular, 


9^. A iUigU{)^, = 


i?f aRl 

aRl i?Q -I- a^Ri * ’ 


(A7) 


which corresponds to the changes 

Ro ^ Ro/\T 2 d\, 7?ii?i|r 2 d|, a ^—a/\T 2 d\ 


(AS) 

or, in terms of the modular parameter (the Teichmiiller 
parameter) T 2 d = ® ^: 

T2d -t -l/T2d- 

On the other hand, 

9,u^ iU29Ul)p^ 

_ f Rl +(a-l)^Rl -ia-l)Rl 
-ia-l)Rl Rl 

which corresponds to the change 
a —t a — 1 (while Rq and i?i are unchanged). (All) 


(A9) 


(AlO) 


The two transformations Ui and U 2 are exactly S and 
T~^ transformations that generate SL(2, Z) (usually used 
in the 2 d conformal field theory literatures), respectively. 

The Euclidean action for the Dirac fermion on this two 
torus is given by 

^E = ^J d^e (det e) i, {T^eYde,) (A12) 

where dete = Y 9 = RqRi, dg^ = and the gamma 
matrices P'^ satisfy {r"^,r^} = 25^^. In terms of the 
space-time coordinates r = RqO^ and x = RY'. 

p27r p27r 

2ttSe = de° d9^ 

Jo Jo 

X ^ (RiYdgo + aRiYdgi + RoYdgi) Ip 

p27tRq p27TRi 

= / dr dx 

Jo Jo 

xij(r°dr + a^r°d, + Yd, 


'll;. (A13) 
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The partition function can be evaluated by the path in¬ 
tegral on the (general) two torus Z{g) = J 
or by the operator formalism 


Z(g) = Tr 




(AM) 


where H' is the ’’boosted” Hamiltonian (in the presence 
of non-vanishing a) corresponding to Se- 

H' = H-ia^P,, (A15) 

Ho 

with 

I n 27 TRi 

H = — J dx'ipT'^d^'tp, 

1 p27TRi 

-Re = ^ y dx (A16) 

being the Hamiltonian and momentum on a ’’flat two 
torus” {a = 0 ). 

The modular invariance for the partition function of 
nonchiral fermions is achieved by summing twisted parti¬ 
tion functions over the spin structures. We thus consider 
the partition function 




E Tra. 

(G^.G;,)GSG 2 




(A17) 


where SG = {1,(—1)-^} is the symmetry group of the 
free fermion theory. Then, the total partition function 
satisfies Z*°\LgL'^) = Z^°\g) for L G SL{2,Z).^^ 


Appendix B: Regularization of the ground 
state-energy 

In this appendix, we regularize the ground-state en¬ 
ergy, which is given by the divergent sum 

Egs[.]{9) = - E (B 1 ) 

sGZ2 + (ax,ay) 


where |s| = y g 2 ^iSj. 

Following Appendix C in Ref. 58, for arbitrary positive 
integer d, we have 


:Z<i + Q, '' ' •> '■>'1 s 

(B2) 


where s,a G |s| = yg^SiSj, ga = det{gdij), and 

Cd^i = - ^ ■ Now we use the equality 

IU; 

= {2Tr)^ - y + 27rn). (B3) 

62 “^ nGZ'^ 


Substituting the above equality, with removing the n = 0 
term in the sum, into (B 2 ), we obtain the regularized 
sum 


H isi'" 


— Cd+ly/^d 

— Cd+ly/ Qd 



E 

n^OeZ'^ 


|y|d+ 




n^OGZ'i 
— 27r2a-n 


27rn 


\d-\-l 


y -h 27rn)e*“ (’'-y) 
(B4) 


Then our regularized ground-state energy is given by 


x=0 


RGSH(ff) = - E 1^1®' 

seZ 2 + (a3;,ay) 

_^ ^—27ria^n. 

= - 03^92 E .3 

n^^OGZ^ |27rn| 


1 ! ; T cos (27ria*ni) 

E —^(B5) 


n^ 0 GZ 2 


(gl^rnrij^ 


Appendix C: Derivation of the claim (36) 

In this Appendix, we confirm the claim (36) by explic¬ 
itly checking how Z^aj(g) transforms under the two gen¬ 
erators Ui = U[M and U 2 of SL(3,Z), defined in Eqs. 
(9) and (13). 

The behavior of Z[a] ( 5 ) under U 2 and C/{ can be directly 
deduced by the properties of the massive theta function 
listed in (5). 

c. Transformation under U 2 : Under 1 / 2 , the metric 

is transformed as in ( 12 ), while the fluxes are transformed 
as (ot, axi tty) —)■ (ot- -I- ay). From (35), we have 

^[G2a](R25R2^) 

= tt ^[ax+/3sy,a^+a^ + (j+i3)sy] (T2d ~ 1"l2Sy) 

Sy^Z-\-(2y 

~ tt ^laj:+/3sy,a^+jSy] (T2d! fl2Sy) 

Sy^Zl-j-Cly 

= ^[a](ff) (Cl) 

d. Transformation under f7(; Under {/(, the metric 
is transformed as in (15), while the fluxes are transformed 
as {ar,ax,ay) —>■ {—ax,ar,ay). From (35), we have 

%(a](R{5C/f) 

~ bl ^la^+jSy, — a^—f!sy]{~^/T2dj1"l2Sy\T2d\) 

Sy^’^-\-CLy 

— bl ®[—Oi—/3sy, —a^— 7 Sy] ('^ 2 dj'Cl 2 Sy) 

Sy^'^-^-dy 

bl ®[ax+/3s„ , 0 ^+ 78 ,,] ('^ 2 dj'Cl 2 Sj,) 

Sy^'^-\-CLy 

= ^[a](ff) 


rrt' 


(C2) 
















13 


e. Transformation under M: Transformation for 
the parameters {Ri, a, /?, 7 } under M is not as obvious as 
the cases of U 2 and U[. We observe that, since the trans¬ 
formation M only involves the change in the x-y plane, 
under M the x- and y- components of the dreibein 
and the metric g^i, (and their inverses) transform as: 

—)■ Mi^e^k, Ca e\ {M ^)k , (C3) 

g., ^ ( 52 )*^' ^ ( 32 )'=', 

where = 1,2, and ( 32 )*^ is defined in Eq. (28). 

To see the behavior of Eq. (34) under M, we first 
note the regularized ground state energy (31) satisfies 
EGS[Ma]{MgM'^) = EGS[al(ff)- On the other hand, the 
second line in Eq. (34) can be expressed as {i,j = 1,2) 

^—27rRQs{s)-\-27riaSj;-\-27ri{oL^-\-'y)sy-\-27riaT 

_ ^—2vRo\/ 32^SiSj+27rz/JoeJ’si-|-27riaT (C4) 


where |u±(s)) are eigenvectors of 

H'(s) = |^a 3 + |^ai+af^ (D4) 

itl 1x2 Kq 

with eigenvalues ±e(s) -I- aSxjRo, where e(s) = 

\J{sx/Rif + (sy/i? 2 )^- Because of P symmetry, 

a^TL'{s)U= TL'is), cr 3 |u±(s)) are also eigenvectors of 
R'is) with eigenvalues ±e(s)-|-aSa;/i?o, and therefore the 
off-diagonal matrix elements are zero, (M+(s)|(T 3 |rt_(s)) = 
(M_(s)|cr3|u+(s)) = 0. 

The diagonal elements, and hence, the transformation 
properties of X±(s) under parity, depend on a choice of 
eigen functions ^±( 5 ). For Sy 7 ^ 0, the following choice 
for the eigenvectors: 


|m±(s)) 


1 

y/2e{s) [£(s) ± Sx/Rl] 


Sx/Ri ± e(s) 
Sy / R 2 


(D5) 


leads to 


where ej* = {a/R q, {a(3+j)/R q)'^. From this expression, 
we can see that the mode-product term in Eq. (34) is also 
invariant under {g,a} —t {MgM"’",Ma}. Therefore, we 
have shown 


Z[Ma](MgM^) = Z[3](5). (C5) 

From the above discussion, we thus confirm our claim 
(36). 


Appendix D: Parity twisted partition functions of 
the surface theory of crystalline topological 
superconductors 


In this Appendix, we explicitly calculate the partition 
functions twisted by parity, which is defined by 

3^ip{x,y)^~'^ = a3ip{x,-y), = 1, (Dl) 

where if is the two-component Dirac fermion. (Remem¬ 
ber that we have doubled the degree of freedom of the 
original theory of Majorana fermions.) Here we define 
y —>■ —y instead y —>■ 2 ttR 2 — y (defined in the main text) 
by parity is just for convenience (the result does not de¬ 
pend on the choice). As mentioned in the text, the par¬ 
ity invariance = H' forces strictly /3 = 7 = 0. 

Then, P acts on the Fourier components of the original 
fermion operators as 

= as-ipis), (D2) 


where s = (s^,, —Sy). On the other hand, the P action on 
the eigen basis x±j defined in Eq. (29), is deduced as 




(m+ (s) 1 0 - 31 W+ (s) ) {u+ (s) 10-31M- (s) ) 
(M_(s)|cr3|M+(s)) (M_(s)|cr3|u_(s)) 


X(s). 


(D3) 


(u+(s)|cr 3 |u+(s)) = (U_(s)|cr 3 |w_(s)) = 1. (D6) 

Alternatively, a different gauge choice 


|m±(s)) = 


1 


V^ 2 e(s) [e(s) ± Sx/Ri 


Sy IE 2 

-Sx/Ri ± e(s) 


(D7) 


leads to 


(u+(s)|cr 3 |u+(s)) = (u_(s)|cr 3 |w-(s)) = -1. (D8) 


In either choice, the result can be summarized as 


’ X+(s) ' 

^-1 ^ 

?7+X-s(s) 

. . 


. V-X-{s) _ 


Sy ^ 0, 


(D9) 


where r]± is an s-independent sign factor. Note that the 
condition = 1 forces 77 ^ = 1. While ri± depends on 
the choice of eigenfunctions, the final results (such as the 
evaluation of the partition functions) do not depend on 
such ambiguity. 

On the parity-invariant line Sy = 0, which exists if 
Oy G Z, the Hamiltonian TLfsxjSy = 0) = ^ 0-3 -|- 
has a ’’chiral decomposition”: 


\ur{Sx)) 


0 ’ 


\ul{Sx)) 


0 


OiR^L € K, 
(DIO) 


which corresponds to the ” chiral eigen basis” xr,l ■ Since 
(Mfl(sa;)|cr3|Mij(s2;)) = -{ul{sx)\'J 3 \ul{sx)) = 1 and 
{uR{sx)\'J3\uLisx)) = {uLisx)\cr3\uRisx)) = 0 (indepen¬ 
dent of the choice of the phases chr/i^), we have 


Xr{sx) 

XlIsx) 




Xr{Sx 

XL^Sx) J ’ 


Sj, = 0, (Dll) 


which does not depend on the normalizations of 
\uRfLisx))- We observe, on the P-invariant line Sy = 0, 
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parity acts like the ’’spin parity” (— 1 )^^, where Fl can 
be thought as the total number of xl[sx) (at Sy = 0 ). 
Thus, we expect that the modular properties of this 
surface theory, as determined solely by the 2 d massless 
modes {sx, Sy = 0 ), will be similar to the modular prop¬ 
erties of the edge theory of ( 2 -|-l)d topological supercon¬ 
ductors protected by Z 2 x Z 2 symmetries. 

/. P-twisted partition functions in the r-direction 
First we evaluate the partition function twisted by P in 
the r-direction, which can be written as 

^^ay 

= Tr 20 .^ r^g27i-i(o^-l/2)Fg-2ii-J?off'l 
^^-2.RoEosYlwf^^iSx), (D12) 


where Oy = 0 , 1 / 2 , Eqs = ^[r](®a;) can be 

written in a pairwise fashion (with respect to P symme¬ 
try): 


In summary, 


^,ay 

const. X Af^ i(T' 2 d)^(" n(T2d) 

X n 0[ iSa-r] (2T2d; ri 2 Sy) 

Sy&Z + 

for Qy = 0 (PBC in the j/-direction), 


= 


const. X 0 [a,„, 2 a,] (2r2d;?’l2Sy) 

s„GZ+ —i 

for Oy = 1/2 (APBC in the y-direction). 

Here the constant prefactors are related to the P eigen¬ 
values of the ground states. 

g. P-twisted partition functions in the x-direction 
Now let us consider the partition function twisted by P 
in the cc-direction. We start with the twisted boundary 
conditions in the x- and ^-directions: 


^[a]('® 2 :) — -^[a]('®a;) X (Sj,) X (Sx), (D13) 

where 

^[a] (S:r) = ^^a^ay ^ CXp | 

+ 2Tri(ar - XK(Sa;)XF(S:r) + xi{sx)XLisx) 

— 2TrSx — 

Ri L 


2'KiaSx 


Xr{Sx)Xr{Sx) - xi{Sx)XL{Sx) 
X^r{Sx)Xr{Sx) + x\{Sx)Xl{Sx)]^ } 


(D14) 


and 


'(/'(a: + 27ri?i, y) = y) ^ 

= e^™"cr 3 V'(a:, -y), 

V'(x,2/ + 27ri?2) = (i^f^)V’(a;,?/)(i^f“=)”' 

= e2”“»i/’(a;,2/). (D18) 

With the above twisted boundary condition, the Fourier 
expansion of the fermion helds can be expressed as 

= E E 

[w+(s)x-e(s) + w-(s)x-(s)] 

+ {2d massless modes}, (D19) 


= n ^exp| 


Sv>0 


X±(s)x±(s) + xi(s)x±(s) 


- 27re(s)^ [x±(s)x±(s) + x±(s)x±(s) 

ill L 


± 


± 27ri^aT- - 

ill L 

2TTiaSx x±(s)x±(s)+X±(s)x±(s) }■ (D15) 


Note that the 2d massless modes {sy = 0) X^j(sa;) would 
be present if Oy = 0. With such pairwise decomposition, 
the 2d massive part for fixed Sy 7 ^ 0 in Eq. (D12) is 
evaluated as 


g27rFo '^(=) n N _ 27ri?o^(s)+27r'iaSx+27riaT- 

Sx 

= 0K,2a,] {T2d] ri2Sy) , (D16) 

while the 2 d massless part is evaluated as 

^fa..a.](^2d)4^,,^_i](r2d). (D17) 


with 

e^”^“X±(s) = ?y±e^”“"X±(s), Sx S | + 0 ,^, = 1, 

(D20) 

where x±(s) are eigen basis of 'H'(s), ^±( 5 ) are the cor¬ 
responding eigenvectors [take the form of (D5) or (D7), 
up to normalization factors], and the term ”2d massless 
modes” is present if Oy G Z. The 2d massless modes are 
given by the sum of the two terms 

e“^Mfi(s^)Xfi(s{(), 

s^GZ-l-aj, 

Y e“’TuL(s{')xL(s^), (D21) 

eZ+Oa; — } 

where Xr,l{^) are eigenbasis of Hfsx^Sy = 0) and 
ur,l{s) are the corresponding eigenvectors in Eq. (DIO). 

From the condition (D20), which relates eigen modes 
with s and s, we only need to take ’’half’ of the degree of 
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freedoms, either modes with Sy > 0 or with Sy < 0 , when 
we calculate the trace for the partition functions. The 
result does not depend on which region for Sy we choose. 
From the above discussion, the 2d massive part for hxed 
Sy ^ 0 in the trace Tr^^ 2 a:„ 
evaluated as 

2-KRoY.s^ez/2+a^ e(s) 




is 


X n 

Sj,GZ/2+ 

~ Q[ 2 aa.,a,-] (”^ 2^/25 ri 2 Sy) 


_ ^—2-n'Ro€{s)+27riaSx-\-27ria^ 


(D22) 


while the 2 d massless part (if present) is evaluated as 

(D23) 


In summary, 
Z 


= 


,ay 

' const. X 

X 0[2a,,a^-a,,] (''■2d/2 + 1/2; 2ri2Sy) 

Sy& + 

for tty = 0 (PBC in the j/-direction), 


const. X n ^[2a^,a^-a,:] (j2dl‘^ + 1/2; 2ri2Sy) 

Sa6Z+-i 

for Ry = 1/2 (APBC in the y-direction). 


Appendix E: Massive modes 01 ™ Ar 2 d',ri 2 ) under 

[flx i^tJ ^ ' 

SL(2, Z) transformations 


In summary, 
Z, 


= 


' const. X 

^ IT Q[2a^ ,aT-] ('^2rf/2; 2ri2Sy) , 

SaGZ+ 

for Gy = 0 (PBC in the y-direction) 


const. X n ©[2a,,,ad {T2d/‘2;2ri2Sy) 

Sj,GZ+-i 

for Oy = 1/2 (APBC in the y-direction). 

(D24) 


h. P-twisted partition functions in the t- and x- 
directions Finally, we calculate the partition func¬ 
tion twisted by P both in the t- and x-directions, 
^ . Using the result from the last section, 

/ ’ f ' ^ 

we now just need to include the additional insertion of 
the parity operator inside the trace. This can be done by 
observing that 


>^X±(s)^ ^=y±X±(s) = e' 
for the massive modes (sy 7 ^ 0 ) and 

'\/ J-. ( o 


2T:iSj- ^ —27rma 


X±(s) (D25) 


XRiSx) 

cai-i _ „ 

Xr{sx) 

_ XL^Sx) _ 

— (73 

. XL^Sx) _ 


for the massless modes (where Sy = 0 as usual). Then, 
the 2 d massive part for hxed Sy 7 ^ 0 in the trace is eval¬ 
uated as 


J27tRq Xls, £ 2/2 + 0, 


— 27ri?.oe'(s)+27r2(a+l)sa;+27r2(aT- — ) 


X Yl l-e 

Sa;G^/2-[-tlx 

~ ^[2aj;,a-r—ax] {'^2d!‘^ “t” 1/2, 2ri25y) , 

while the 2 d massless part is evaluated as 

^K,ad(^ 2 d)A^^_l_^^_l](T 2 d). 


(D26) 


(D27) 


In this Appendix, we discuss how the (products of) 
massive modes ©J~™ ](T' 2 dUi 2 )) dehned in Eq. (65), 
transform under SL(2,Z) generated by U[ and U 2 - This 
can be deduced from the modular properties (5) of the 
massive theta functions with modular parameters T 2 d, 
2 T 2 d, T' 2 d/ 2 , and T 2 d /2 -1-1/2 (we denote the mass param¬ 
eter m = ri 2 Sy in the following equations): 

(i) For 0[a,,a,](T-2d;™): 

0[a,,a,l(r2d;TO) ^ 0[a,.a,] (“ l/r2d; TO|T2d |) 

©[—a,,a,] i.^2d^ ^) i 

Q[a^,a^]{T 2 d;m) ©[a, .a,] ('r2d + 1 [ w) 

= ©[a„a,+a,] (T2d; m) ; (El) 

(ii) Eor 0[a^_a^](2r2d;m): 

©[a,,a,]( 2 T 2 d; W) -% ©[a,,a,] (-2/r2d; w|T2d| ) 

= ©[-a,,a,] (T 2 d/ 2 ; 2 TO) , 

©[a,,a,] (2T2d; m) ^ ©[a,,a,] (2T2d + 2; ITl) 

= ©[a,,2a,+a,] (2T2d; w) ; (E2) 

(hi) Eor 0[a^_a^](T2d/2;2TO): 

©[a,,a,] {T2d/2;2m) ©[o,.a,] (“ l/2r2d; 2TO|r2d |) 

= ©[-a,,a,] ( 2 T 2 d;m), 

©[a,.a,l (T 2 d/ 2 ; 2m) ^ ©[a,.a,] (T 2 d /2 -f 1/2; 2m); 

(E3) 

(iv) Eor 0[a^_a^](T2d/2-|-1/2;2 to): 

©[a,,a,] (T 2 d /2 -f 1/2; 2m) 

©[a,,a,] (-l/ 2 r 2 d -f 1/2; 2 m|r 2 d|) 

©[—a, —2a, ,a,+a,] ('^ 2 d /2 1/2, 2m) , 

©[a,,a,] (T 2 d /2 -f 1/2; 2m) ^ ©[a,,a,] {T 2 d /2 + 1; 2m) 

= ©[a,,a,+a,] {T2d/2; 2m). (E4) 
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Therefore, 


[ax 



[ a-T-, (lx ] ^ 


^04 

^[dx,dx+dx]’ 

[ax 



r\iii 

[ a-T-, (lx ] ^ 


U~^ 

^ ^[dx,2dx+ax] ■ 

Ail 

'k 

,“t] 

ul 

[ a-T-, dx ] ^ 

r\iii 

f/-i 

^2 p.tv 

[dx idx] ^ 

Av 

[ax 

,aT] 


piiv 

'^[—dx—2ax 

,a^+a^] I 

U~^ 

'-'[dx,axl ^ '^[dx 


(E5) 


The constant prefactors are again related to the P eigen¬ 
values of the ground states, which can be absorbed to the 
(redefined) weights as we consider the partition sum. 


The total partition function is then given by 


Appendix F: SL(2,Z) invariance of the total partition 
fnnction for % = 1/2 


The parity-twisted partition functions for ay = 1/2, as 
computed in Appendix D, are summarized as follows: 

a„ = i = const. X 0[a,„,2a^]('^2d; ^' 12 ), 

= i = const. X 0[2a^,a^]('^2d; ^' 12 ), 

= i = const. X 0[2a„,o^-a^]('’’2d;’’l2), 

(FI) 

where we have introduced ^{T 2 d',rl 2 ) as: 

ydx ,o,T \ ' 

®[a.,a.]= n [ 0 [a..ad(^ 2 d;ri 2 S,,)]", 

s„ez+-i/2 

®[L,a.] = n 0 [a..ad( 2 'r 2 d;ri 2 Sy), 

s„GZ+-l/2 

0 [r,,a,] = H 0 K.ad('C 2 d/ 2 ; 2 ri 2 Sj;), 

s„GZ+-l/2 

0[L.a.]= n 0K.ad(W2 + l/2;2ri2S,). (F2) 

s„GZ+-l/2 


= C10[O.O] + ^20[O,1] + C30[1,O] + ^40[1,1] 

+ 2(e50[Q 0] -I- e60[i_o] + C70[o)o] + C80Jo*i] 
+ C90[O,O] + ^10®[0,i]) ■ 


From the modular properties of 0 (and thus of 0* ™) 
discussed in Appendix F, we can see that Zf°*_i,(gp) 

can be made SL(2,Z) (generated by U[ and U 2 ) invariant 
for any number of Dirac fermion flavors, N, if we choose 
Ci = 1 for all i (more precisely, we just need £2 = £3 = £4 
and £5 = ... = £ 10 ). 
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